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Abstract 

We analyse and numerically simulate the full many-body quantum dynamics of a 
spin-1 condensate in the single spatial mode approximation. Initially, the condensate 
is in a "ferromagnetic" state with all spins aligned along the y axis and the magnetic 
field pointing along the z axis. In the course of evolution the spinor condensate 
undergoes a characteristic change of symmetry, which in a real experiment could be 
a signature of spin-mixing many-body interactions. The results of our simulations 
are conveniently visualised within the picture of irreducible tensor operators. 
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1 Introduction 



This paper is a tribute to Bruce Shore's role in establishing the concept of 
atomic coherence and the related techniques of irreducible tensor operators 0. 
We demonstate that this viewpoint emerges quite naturally when describing 
the evolution of a spinor condensate. During the past decade condensed Bose 
gases have energed as flexible test system to explore the rich structure of 
many-body physics. A systematic introduction into cold gases in general and 
spinor condensates in particular may be found, e.g., in the book by Pethick 
and Smith . With the development of purely optical traps jli 0, , we are 
in the position to hold atoms with their complete hyperfine submanifold, thus 
entering the regime of spinor physics f^H,0,H,|^. This is exactly the condition 
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under which one expects the techniques of atomic coherence [l|, E3, lUl to be 
most helpful. 
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There are two major approximations used in the hterature when describing 
spinor condensates theoretically: the common mean-field approximation (MF) 
and the single mode aproximation (SMA). Within the SMA, the spatial depen- 
dence of the Bose field operator describing the spinor gas is assumed given and 
decoupled from the internal dynamics, as illustrated schematically in Fig. 1. 
The reason one can separate the local redistribution of particles amongst hy- 
perfine sublevels due to collisions from the long range dynamics is clear: a 
typical trapped alkali atom is paramagnetic, but provides only a single Bohr 
magneton of magnetic moment. Consequently the resulting magnetic dipolar 
interaction is very weak compared to the typical interactomic Van-der-Waals 
energy. 




^(x) = <l>{x) (g) {a_i, ag, a+l) 

Figure 1. Schematic diagram of the harmonically trapped spin-1 ^^Rb condensate 
in the single-mode approximation. 



The single mode approximation is commonly introduced on top of the mean- 
field approximation. One exception is the paper by Koashi and Ueda 2| where 



the SMA is introduced directly into the many-body treatment (see also 
This, and also neglecting the quadratic Zeeman part of the energy, allowed the 
authors to construct an explicit solution to the problem including collisions. In 
this paper we include the quadratic Zeeman part which plays a crucial role in 
the dynamics of the condensate as realized with '^^Rb or ^^Na atoms [1,0! M 



Consequently, we are left with an interacting three mode Hamiltonian, whose 
dynamics will be the subject of study in this article. 

It should be pointed out that our approach would not apply directly to the 
recently realised Chromium Bose-Einstein condensate [l5| where the dipolar 
interaction is comparable to the Van-der-Waals energy. Even less so in the case 
of heteronuclear molecules pursued by many experimental groups around the 
world where the dipolar interaction is expected to be the dominant mechanism. 
Independence of the internal and spatial motion is completely lifted e.g. in the 
proposed Einstein-de Haas effect [3 ■ The SMA is obviously not aplicable 
under the conditions for condensate fragmentation and domain formation (see 
[lil and references therein). 

The article is organized as follows: In Section 2 we review the derivation of the 
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three mode spinor Hamiltonian. For later reference, we will also discuss the 
standard mean-field approximation in the language of Hamiltonian mechanics 
in Section 3. Visualizing the results of a full many-body evolution of the three 
mode system is a nontrivial problem. Therefore, we introduce the concepts of 
irreducible tensor operators in Section 4. Finally, we discuss and compare the 
results of the many-body calculation in Section 5. 



2 Many-body approach to spinor condensate dynamics 



The F = 1 bosons are described in the usual way by a field operator with the 
three spherical components, ipkif), where k = 1,0,-1 (in this order) labels 
the z projection of the atomic spin. The effective low energy Hamiltonian of 
such bosons in a homogeneous magnetic field is written as: (cf. 0, H, EHI3|) 

Here, 

^0 = E / d'r^plir) [ - + [/(r)]4(r') (2) 

k 

contains the kinetic energy and the trapping potential which holds the atoms 
in place, and 

nB = T.Ek [ d'ri,i{r)Mr) = E ^kA (3) 
k •' k 

describes the effects of the magnetic field. For weaker magnetic fields, this is 
mainly the linear Zeeman effect, Ei — Eq Eq — E^i oc 5, and the quadratic 
Zeeman, or Paschen-Back, effect, Ei + E^i — 2Eq oc B^. The latter is crucial 
for the purposes of our analyses. The general Breit-Rabi formulas for Ek may 
be found in Ref. l2C 



Both Ho and Hb describe single-body physics. Two-body collisions enter via 
H^cli S'^d ^coii- The former describes the famihar spin-independent density- 
density interaction, 

^Si = |E/ d'ri.l{r)^j{r)Mr)Mr) = |/ rfV:n^(r):, (4) 

k,l 

where n(r) = J2k^k{r) = Z^fc '^fc(^)V^fc(^'), and :: stands for the normal order- 
ing of the field operators. The latter describes the spin-dependent part of the 
collision, 

n'l='jjd^r:jr\r):, (5) 
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where ^{r) is the density of the angular momentum of the condensate, 



kl 



(6) 



The Cartesian components of the single-particle spin operator, F, are defined 



as: 
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0-1 



(7) 



The coUisional interaction constants, Cq and C2, are expressed in terms of the 
scattering lengths and 02, for two spin-1 atoms in the combined symmetric 
channels of total spins, respectively, and 2, as \^ 

_ ATih^{aQ + 202) _ A'nfi^{a2 - ao) 

^0 — 7rv7 5 ^2 



3M 



3M 



A typical value of C2 is one-two orders of magnitude below cq. This leads to a 
natural separation of time scales: on smaller time scales, one may neglect Ticoii? 
whereas for longer time scales it may be of crucial importance. The ultimate 
example of a longer time scale is the ground state of the condensate, which 
may be polar or ferromagnetic depending on the sign of C2 An example 
of the opposite are experiments where the spinor condensate is prepared in a 
^articular magnetic state, and the evolution of the condensate is then observed 



particular 



|. In this case, an approximation commonly used in the literature is 
the so-called single mode approximation (SMA). Technically, SMA is recovered 
by assuming that the spatial dependence of the field operator is fixed, ipkij") = 
(j){r)ak, J (i^r|0(r)p = 1. Under the SMA, the Hamiltonian of the system of 
spin-one particles becomes: (cf. ET 



n 



-p:f. + qP2 + 3-1 ( 



2Af 



(9) 



Here, 92 = 02! ciV|0(r) |^ p = - E,)/2, and g = (^1 + ^-i)/2 - Eo. The 
"wide hat" symbol denotes mapping of a single-body operator into the Fock 
space. 



k,l 



(10) 



with ai being the mode annihilation operators; in our case, / = 1, 0, — 1 (in this 
order). The "wide hat" should not be confused with the "normal" hat which 
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just signifies that a quantity is an operator in the Fock space. In particular, 
= is the full angular momentum operator, and A/" = X3 is the Fock-space 
operators of the total particle number, where X3 is a 3x3 unit matrix. The 
second line in (9) is found using that, 

:i^':=^'-P = i^'-2A/'. (11) 



In numerical terms, constituents of (9) are sparse matrices; these are handled 
very efficiently by common high-level numerical packages. Some details on how 
these were implemented numerically may be found in Appendix A. 



3 The canonical formulation of the mean-field dynamics 



3. 1 Introducing the mean field 



In the mean field approximation, operators are replaced by complex num- 
bers Cfc- However, this transition is not uniquely defined. Using different or- 
derings (normal, symmetric, etc.) of creation and annihilation operator yields 
different functional forms of the quantum Hamiltonian. In quantum mechan- 
ics these functional forms represent the same Hamiltonian. In the mean field 
picture where ordering is immaterial for the complex numbers Cfc, Cki these dif- 
ferent functional forms result in different classical Hamiltonians. In the scalar 
case, the corresponding corrections are of the order of 1/A^, where N is the 
total number of atoms. For the spinor condensate, these corrections are of the 
order of one over the population of Zeeman sublevels. There can well be patho- 
logical cases when one level is much depleted, leading to large corrections. It 
is therefore of importance to have the procedure of introducing the mean field 
approximation for the spinor condensate clearly defined. 

There is a number of ways the mean field picture can be introduced. The 
most rigorous one is via the phase-space technique (see, e.g., jl^l and refer- 
ences therein), where the time-dependent Gross-Pitaevskii equation appears 
on dropping quantum noises in the equations for the phase-space amplitudes. 
This procedure is equivalent to using the normally-ordered form of the inter- 
action Hamiltonian, replacing the field operators by the c-numbers, and then 
postulating the resulting c-number function to be the classical Hamiltonian. 
The lack of rigour is clearly seen in the fact that this procedure leaves it un- 
clear how quantum averages are expressed by the c-number amplitudes. Some 
insight can be gained from the fact that if we postulate the many body state as 



a coherent state, 
the condition 



C) = ICi) ICo) IC-i)) and then minimize the energy subject to 
J\f C) = IICIP = ^) the time-independent Gross-Pitaevskii 
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equation is recovered. This recipe is in fact general: a correct mean-field ex- 
pression for the average of an operator X is found as (^x'^ ~ (*» This is 
equivalent to replacing field operators by c-numbers in the normally ordered 
representation of X. 



3.2 The classical Hamilton equations of motion 



Following the above recipe, we postulate the classical Hamilton function to 
be, 



-p{\C.\'-\W)+<i{\Ci? + \C-i\') 

92 



+ 



ICi 



ic-ir-2icinc-ir+2iconcii' 



+ 2|Conc-ir + 2CocrCi + saxiC-i 



(12) 



where (9) was used to find the explicit form of the Hamilton function. To use 
the powerfull methods of analytical mechanics we have to know the canon- 
ical conjugated momentum tt to the coordinate C- If) by definition, x and 
p are three-dimensional canonical conjugated variables, defining the complex 
variable C 

x-F ip 



V2 



(13) 



than we can identify tt = i^* . One can verify easily, that this is a proper 
canonical transform, as 

with the following conventional definition of the Poisson brackets. 



So we obtain the Hamilton equations of motion 

^0-{0,i^(C,7r)}^,,,^-^-— , 



(15) 



(16) 



dt 



dCi 



By employing complex coordinates we have doubled the dimension of the phase 
space making the two sets of equations of motion redundant. Explicitly, we 
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find the three components of the time-dependent Gross-Pitaevskii equation 



Ho 



-p + qO 

p + q 



Hi 



C-iCo 





ICiP + IC~iP 

crco 





CiCo 



(17) 
(18) 



with = iCip- 
Sengstock group 



^, the 



^ For the initial condition as in the experiment of the 
these nonhnear equations have been solved analytically 



in terms of the periodic Jacobi elliptic functions, see 



3.3 Comparing the many-body and mean-field dynamics 



Consider now how an effective mean-field amplitude is recovered in the many- 
body techniques. Were the population and coherence properties of the matter 
field investigated in a measurement, all the relevant information enters through 
the single-body density matrix. 



Pik{t) 



dl(t)di(t) 



(19) 



where |\l/o) is the initial many-body state, and di{t) is the Heisenberg field 
operator. 



di{t) = U{ty diU{t), U{t) = exp{-im). 
It can be also written as an object in the Schrodinger picture. 



Pik{t) = hoit) 



^o(t) 



(20) 



(21) 



where |\l/o(t)) = W(t) l^&o)- The order or matrix indices was chosen so as to 
have a natural property. 



^^o(t) A ^o(t)) = Trp(t)A (22) 

for any single-body operator A. The "wide hat" symbol is defined by Eq. (10). 

In the mean-field approach, the density matrix is expressed by the mean-field 
amplitude, Cit): 



(23) 



Mathematically speaking, this relation stipulates that has a single nonzero 
eigenvalue equal to A^, with N~^^'^<^(t) being the corresponding eigenvector. It 
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is then only natural to introduce effective mean-field amplitudes as eigenvec- 
tors of p(t): 



Pikit)= E rs{t)C\'\t)Ci^\t), C^'^-C^''> = NSss'. (24) 

s=l,2,3 



In general, we have three eigenvectors of p(t), N~^^'^(^^'^\t), and three related 
eigenvalues, Nrs(t), ri{t) + r2{t) + rsit) = 1. For a mean- field approximation 
to be recovered, there must exist a dominant eigenvalue, ri{t) ^ r2{t),r^{t). 
Below in Section 5 we shall see that it is exactly the situation occuring at early 
stages of the condensate evolution. At later stages all three eigenvalues are of 
the same order of magnitude so that the mean-field approximation breaks 
down completely. 



4 Time evolution of atomic coherence 



4^.1 Irreducible tensor operators 



The physical information contained in the correlations expressed by the single- 
body density matrix is atomic coherence 0. A natural framework for dis- 
cussing such correlations is the picture of irreducible tensor operators [l[ E^, 



An irreducible tensor operator of angular momentum j is a set of 2j + 1 
square matrices, Tjm, m = —j, — j + 1, ■ ■ ■ , j. Under three dimensional rota- 
tions, they obey the same transfomation laws as spherical harmonics Yjm{0, (p), 
thus forming a natural basis for expressing any physical quantity with angular 
momentum j. Furthermore, irreducible tensor operators form a full orthonor- 
mal set in the space of matrices in respect of the Hilbert-Schmidt norm: 



Tr Tj^T^t^, — Sjj'Smrn'y (25) 



and for any (2 J + 1) x (2 J + 1) matrix A, 

2J j 

^ = E E T,„, cf = Tr T^A. (26) 



=0 m= 
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4-2 Dipole and quadrupole momenta of the spin-1 condensate 



4-2.1 Cyclic basis 

For the spin-1 condensate, the single-body density matrix (19) is a 3 x 3 matrix 
and hence only contains angular momenta 0, 1, and 2: 

p{t) = cgToo + ^TTim + E c^T2m. (27) 

m=— 1 TO=— 2 

Cq and have a very simple interpretation: Cq is in essence the total popula- 
tion, 

and c^ is related to magnetisation, 

= (29) 

where the spherical components of the angular momentum operator are de- 
fined as, — ^ (^^x ± i^y^ /V^, = J^z- C21J. characterise the quadrupole 
momentum (alignment) of the atomic state. 



4-2.2 Cartesian basis 

It is instructive to rewrite p{t) in the Cartesian basis, (with a, (5 — x, y, z) 



Paf^it) = (^0 a}{t)d^{t) ^0/ , ^g^^ 

a± = T {dx =t iay) I \p2., do = dz- 

Then, atomic population corresponds to the trace of the density matrix, 
magnetisation — to its antisymmetric part, and alignment — to its traceless 
symmetric part. Taking the symmetric part of an Hermitian matrix is the 
same as taking its real part, hence contains a mix of the population and 
alignment while S^p is pure magnetisation. Explicitly, 

Pap = ^PaP + '^PaP = ^Pafi + 9 E ^afiaJ^cr, -^a = E ^oc^aPaP, (31) 

T a,P 



where J- = l^!Fj is the magnetisation vector, and Sapa is the fully antisym- 
metric tensor, Sapa = —spaa — —^aap, ^xyz = 1- Introducing the quadrupole 
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tensor as per Ref. [11 



TV 

Qa/3 = ~ ^PajS, (32) 

we find the irreducible-tensor-operator expansion of p{t) as 

jV i 

Pap = ir^aP + 7J 51 ^apa^a — Qafi- (33) 



4-3 Atomic coherence under the mean-field approximation 

In tlie mean-field picture, = Ca(p, where the Cartesian components of 
C = C + K" found in full analogy to (30). The population-and-alignment 
part of the density matrix is then, 

^Pap = CCp + CQ. (34) 

Magnetisation also follows easily, 

:r=-2[C'xC"]. (35) 

We can assume that either C, is real, = 0, or the real and imaginary parts 
of C, are orthogonal, C,' _L C," ■ Indeed, since we can always find such that 

3fJCe*'^ ■ 53Ce*'^ = C' ■ C" cos 2^ + ^ ^ ^ sin 2^ = 0, (36) 

this orthogonality can be achieved by an overall phase transformation. Then, 
Eq. (34) explicitly diagonalises 3ftp, making it evident that in the mean-field 
approximation 3fJp always has a zero eigenvalue. Another two eigenvalues equal 
C,'^ and C"^. If C, is real, two eigenvalues of are zeros, leaving the system in 
a maximal-alignment state along In this case, the magnetisation which is 
related to '^p vanishes. The inverse is also true: zero magnetisation is charac- 
teristic of real C,. 



An immediate word of caution is necessary here. Orthogonality of 3?<^ and 
is not preserved by the Gross-Pitaevskii equation (17). One can rewrite (17) 
so as to maintain the said orthogonality, but the resulting equations are rather 
bulky and do not seem to allow for a new insight; worse still, under certain 
conditions the phase evolution is discontinuous. 
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4-4 Equations of motion for the atomic coherence 



In order to gain further insight into the evolution of the atomic coherence, 
consider the equations of motion for the single-body density matrix written 
in terms of the magnetisation and quadrupole matrix. Direct calculation with 
Hamiltonian (9) yields, 



a 

QaP = (pQ<T(3 + ^SfSzJ^a ) + {« ^ /?} (37) 



+ 



Q 



a/3 



coll. 



The benefits of the representation of irreducible tensor operators are now 
clearly seen. Collisions are spherically symmetric so that only the linear and 
quadratic Zeeman effects may couple magnetisation and alignment. In actual- 
ity, additional reflection symmetries leave only quadratic Zeeman effect to do 
it. The linear Zeeman effect leads only to the overall Larmor precession around 
the z axis. Furthermore, collisions only contribute directly to the evolution of 
the quadrupole. 



Q 



coll. 



(3J 



In general, this contribution contains nontrivial two-body correlations and 
cannot be simplified any further. Within the mean-field approximation. 



Q 



al3 



coll. 



^92 



(39) 



Note that here we do not assume that ^' ± 



Equations (37) allow one to make qualitative statements about the evolution of 
the system's coherence independently of the details of collisional interactions. 
Assume that the initial state of the system is as in experiment j^. Then, the 
initial magnetization vector is along the x axis, JF^ = 0. Since JF^ = 0, cf. Eqs. 
(37), magnetization ^ never leaves the xy plane. As a result, under the mean 
field approximation the alignment eigenvectors stay in a plane orthogonal to 
the xy plane, cf. (35). If at a particular time magnetisation vanishes, the mean 
field picture predicts that the system must be axially symmetric. Hence, under 
the SMA, absence of axial symmetry at zero magnetisation is a clear signature 
of many-body effects. Since thermalisation via symmetric collisions can only 
increase the symmetry of the system, it appears plausible that this statement 
should survive beyond the SMA. Verifying this conjecture is subject to further 
analysis. 
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5 Results and discussion 



5. 1 Fock space of the spinor condensate 



The choice of the basis in the Fock space is inherently linked to the angular- 
momentum properties of Hamiltonian (9). A native bookkeeping is to intro- 
duce the basis states by the populations of the single-particle states with 
Fz = 1, 0, —1. By construction, such basis state in the Fock space is an eigen- 
state of the three mode population operators: 

A \Nu No, N^i) = Nk \Nu No, N.^) , = 1, 0, -1. (40) 
These states also happen to be the eigenstates of A/" and JF^,: 



Af \N,, No, iV_i) = N\N,, No, N_,) , 
|iVi, A^o, AT.i) = M \Nr, No, N.^) 



(41) 



where A^ = A^i + A'o + A^-i and M = Ni - N_i. Since H commutes with Af 
and J-'z, states with different A^'s ans M's are not mixed by free evolution. 
Another simplification when using this "natural" basis is the ease with which 
the many-body operators are implemented (Appendix A). 

In the absence of the quadratic Zeeman effect (QZE), Hamiltonian (9) con- 

- 2 

serves the full angular momentum, ^ , and its component, J-'z- In this case 
the time evolution is solved ^3] by simply classifying the basis states by the 
eigenvalues of the three operators: 

^ I A^, F, M)^ = N \N, F, M)^ , 

J'^\N,F,M)^ = F{F + l)\N,F,M)^, (42) 
:Fz \N, F, M)^ = M \N, F, M)^ . 

Here, < F < A^ and —F < M < F. The subscript distinguishes the 
"angular momentum" basis from the "natural" basis introduced above. 

In the presence of the QZE, F is not a good quantum number any more, 
while M is still conserved. Hamiltonian (9) can then be diagonalised within 
subspaces of given M of the "natural" basis. Conversely, rotating the states 
by the radio-frequency field as in Ref. jll conserves F while mixing different 
M's. Using the "angular momentum" basis may thus be of some help. In the 
calculations summarised in this paper, the natural basis was used. Exploring 
the advantages of the angular momentum basis remains a subject to further 
work. 
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5.2 Numerical procedure and the choice of parameters 



All simulations summarised in this paper were carried out with the experiment 
by Sengstock and co-workers in mind In this experiment, the condensate 
is firstly prepared in a state where all atomic spins are antiparallel to the z 
axis. A radio-frequency pulse is then used to re-orient the spins antiparallel 
to the X axis. In our notation, this state looks as 



where = ti/2. Since the manifold of M = — is non-degenerate, |0, 0, A^) is 
also an eigenstate of the total angular momentum, |0,0,A^) = |A^, A^, — A^)jc-- 
Furthermore, since the evolution does not mix different M's, |0, 0, A^) must 
be an eigenstate of the Hamiltonian, providing us with a test of numerics. 
Another test is setting d = tt, yielding l^'o) = |A^5 0,0) = |A^, A^, A^)^, which 
is also an eingenstate of the Hamiltonian. Having thus checked that both for 
■(9 = and for i) = it "nothing at all happens," we proceeded with a series of 
simulations with (9 = 7t/2. 

To compare the many-body results to the mean-field approximation, we also 
ran a mean-field simulation of the system using (17). It easy to see that, with 
our choice of the many-body initial state, Eqs. (23) are exact at t = 0. This 
provides us with the initial condition for the mean-field simulation. This is 
certainly the same initial mean-field amplitide as used in Refs. [sl. 11^. 

Simulations presented in this paper were run for the following set of parame- 
ters: p = 0, q = 1, and N q-), = —0.3. These differ from the parameters charac- 
teristic of experiments jsi Il4|. p ^ q, Ng2 = — 0.03g, by setting p to zero, and 
enhancing nonlinearity by an order of magnitude. It should be emphasised that 
observation times in these experiments are too low for the many-body effects 
to manifest themselves. Simply extending the "observation" time in simula- 
tions allows for the many body effects to be seen, but produces two many 
oscillations of physical quantities on the relevant time scales. Enhancing the 
nonlinearity allows one to produce comprehensible graphics, while making no 
qualitative difference to the condensate evolution. Setting p to zero is of no 
consequence: in essence this means introducing a reference frame following the 
Larmor precession. 

The number of particles in our simulations was A^ = 50. Importantly, this 
number is large enough to clearly separate the relaxation and the revival time 
scales. Indeed, relaxation is exacty that what we expect from a many body 
treatment, whereas revivals are artifact of the model. At the same time, leaving 
the atom number relatively low simplifies numerics enourmously, making the 
problem easily tractable on an average desktop computer. 




(43) 
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Figure 2. Evolution of the condensate magnetisation as function of the scaled time. 
A'^ = 50, p = 0, q = 1, Ng2 = —0.3. Only the x component of the magnetisa- 
tion is nonzero. In the mean-field approximation, magnetisation oscillates without 
relaxation. 




Figure 3. Scaled eigenvalues of the single-time density matrix. = 50, p = 0, q = 1, 
Ng2 = —0.3. In the mean field approximation, one eigenvalue always equals one and 
the other two eingenvalues are zero. 



5.3 Evolution of atomic coherence 



Figure 2 shows the evolution of the condensate magnetisation, Eq. (35), calcu- 
lated both within the many-body model and under the mean-field approxima- 
tion. With our choice of the inital state and parameters, only the x component 
of the magnetisation is nonzero. In the mean-field approximation, this com- 
ponent oscillates periodically between values —N and A^. Exact periodicity 
follows from the analytical solution 0]. The many body result shows decaying 
oscillations with the same period. Qualitatively, the decay can be understood 
if considering the simplest generic model of collisional nonlinearity: the Kerr 
oscillator [23[. The latter is a single-mode quantum oscillator with quartic 
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nonlinearity, 



HeS = uJesn+^n{h-l). (44) 

Assume the Kerr oscillator is initially in a coherent state with amplitude a. 
A straightforward calculation of the average position of the oscilator yields, 

{x{t)) = -^3f?{aexp [ - iu^^t + iVeff(e-'«««* - l)] } 

«^XMFWe-^««^eVV2^ befTi|«l. (45) 

where A^eff = 

XMFit) = -Ls)^[tte-^('"'^«+^««^<^ff)*] (46) 

is the mean-field approximation to {x{t)) which accounts for the frequency 
shift but not the dephasing. Given the obvious crudeness of the Kerr model, it 

reproduces the evolution of the magnetisation, Fig. 2, with surprising accuracy. 
Setting Ncsg^Q ~ 2A^(y'| results in a good quantitative estimate of the decay 
time, Tdecay ~ {^fN 92)^^ ~ 23. The Gaussian nature of the decay is also 
evident in Figs. 2 and 3. 

Contrary to what could be expected, the system does not evolve towards 
a spherically symmetric state. A formal indication of this may be seen in 
Fig. 3. This figure shows scaled eigenvalues of the single-body density matrix, 
cf. Eq. (24). In the mean field approximation, one eigenvalue always equals 
one and the other two eingenvalues are zero. This holds for the initial state, 
but for nonzero times wc find that all three eigenvalues are nonzero. Most 
interestingly, they do not evolve towards the value of 1/3 indicative of spherical 
symmetry. Two of them stabilise at ^ 0.22, while the third one remains close 
to 0.56. Since the magnetisation decays completely, cf. Fig. 2, the persistent 
asymmetry must be of quadrupole nature. 

Being a mathematically rigorous way of characterising alignment, the quadrupole 
tensor is not the best choice so far as visualisation of the results is concerned. 
A much more convenient way is thinking in terms of the alignment ellipsoid. 
The latter is found by diagonalising the real part of the density matrix, cf. 
Eq. (31). The eigenvectors of are the alignment axes. The eigenvalues of 
are positive, their sum equals population, qi + q2 + (ld, = A^, and differences 
between them characterise alignment. In particular, three equal eigenvalues, 
0.1 — 0.2 — % — N/3, mean the atomic state has no alignment at all. The 
system of three orthogonal eigenvectors and three positive eigenvalues is nat- 
urally visualised as a three-dimensional ellipsoid. If we additionally shift the 
centre of this alignment ellipsoid by the magnetisation vector, the resulting 
picture contains all details of the atomic coherence, cf. Figure 4. An additional 
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Figure 4. Visualisation of the time evolution of the condensate coherence. The pic- 
tures shows the alignment ellipsoid, the centre of which is shifted by the magneti- 
sation vector. One of the alignment axis is always collinear with the magnetisation 
vector, so that the only parameter specifying the ellipsoid's orientation is the tilt 
angle. In the mean-field approximation, the ellisoid reduces to a fiat ellips orthog- 
onal to the magnetisation vector. Top row: initial stages of the evolution at which 
the mean-field picture is valid. Left: the initial state of the condensate coherence. 
The alignment ellipsoid is reduced to a disc. Right: the condensate coherence at zero 
magnetisation. The alignment ellipsoid is reduced to a straght line in the yz plane 
with 7r/4 tilt. Bottom row: later stages of the evolution at which the mean-field 
picture fails. Left: a typical picture of coherence at intermediate stages of the evo- 
lution. Right: the final state of condensate coherence. The magnetisation and tilt 
are zero, and the ellipsoid is axially symmetric in respect of the quantisation axis. 

simplification is that the smallest axis of the eUipsoid was always found to be 
collinear with the magnetisation vector. This also holds if p 7^ when the 
whole picture is subject to the Larmor precession. Thus the only parameter 
we need in order to fully specify the orientation of the alignment ellipsoid is 
the tilt angle; the latter is defined as the angle between the largest axis of the 
ellipsoid and the quantisation (z) axis. Together with the eigenvalues of [cf. 
Eq. (31)] and the magnetisation vector, this angle completes characterisation 
of the atomic coherence. 

Visualisation of atomic coherence at different stages of the evolution may be 
seen in Fig. 4. Evolution of the afignment eigenvalues and of the tilt angle is 
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Figure 5. Evolution of the quadrupole. A'' = 50, p = 0, q = 1, Ng2 = —0.3. Left: the 
alignment eigenvalues, right: the tilt angle. The smallest eigenvalue corresponding 
to the alignment axis along the magnetisation is shown as solid line, the other 
two — as dashed lines. The modulus of magnetisation (dotted line) helps to visually 
synchronise this picture with Fig. 2. In the mean field approximation, the smallest 
eigenvalue is exactly zero while the other two eigenvalues and the tilt oscillate 
periodically without relaxation. 

shown in Fig. 5. At early stages of the evolution, the mean field picture holds 
with a very good accuracy. In the mean field approximation, the alignment 
ellipsoid reduces to a flat cllips orthogonal to the magnetisation vector. There 
are two characteristic moments of the mean-field evolution [cf. Eq. (35) and 
Fig. 5, top row]. If the magnetisation is maximal, the alignment axis are equal; 
the alignment ellipsoid becomes a disk (Fig. 5, top- left). If the magnetisation is 
zero, only one alignment eigenvalue is nonzero; the alignment ellipsoid reduces 
to a straight line at 7r/4 tilt (Fig. 5, top-right). The mean-field evolution of the 
atomic coherence is in essence periodic oscillations between these charateristic 
configurations. 

On later stages of the evolution the mean field picture breaks down. All three 
eigenvalues of are nonzero, so that the alignment ellipsoid acquires thick- 
ness, becoming an ellipsoid proper (Fig. 4, bottom-left). The tilt oscillations 
are damped and finally cease (Fig. 5). The final state of atomic coherence is 
depicted in Fig. 4, bottom-right. The magnetisation is zero so that the eUip- 
soid is centered at the centre of coordinates. The larger axis is along z, the 
other two are equal so that the alignment ellipsoid is axially symmetric. This 
state of atomic symmetry is impossible to obtain in the mean-field picture. 

Thus in the course of the evolution atomic coherence undergoes a complete 
change of symmetry. Initially, the system is prepared in an axially symmet- 
ric state in respect of the At intermediate stages of the evolution, all 
symmetries are broken. Finally, the system again tends towards an axial sym- 
metry. Not unexpectedly, the final symmetry coinsides with the dynamical 
symmetry, with the symmetry axis parallel to the magnetic field. 
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6 Conclusion 



We have modelled the quantum many-body dynamics of a spin-1 condensate 
numerically within a single mode approximation, using the representation of 
irreducible tensor operators so as to conveniently visualise the condensate evo- 
lution. Comparing results of the numerical simulation to those obtained within 
the mean field approximation shows a number of characteristic features of the 
many-body evolution which are impossible to obtain under the mean field ap- 
proximation. In particular, the final state of the condensate coherence is of a 
characteristic quadrupole nature, axially symmetric in respect of the quanti- 
sation axes, whereas this kind of atomic coherence is impossible in principle in 
the mean field approximation. Experimentally observing this state of atomic 
coherence would be a distinct signature of spin-mixing collisions without the 
mean-field approximation. 
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A Implementation of many-body operators in the Fock space. 

To be specific, consider how mode annihilation operators, dk, are implemented 
in a high-level numerical package. Such packages handle sparse rectangu- 
lar matrices very effectively, so the first step is to replace a triple index, 
{iVi, iVo, A^-i}, characterising a basis state in the Fock space, by a linear 
index. A helpful trick is to start from introducing an extended index, n = 
1 + iVi + (A^ + l)A'o, where < A^i, Nq < N. The third index, N_i, is found 
from Ni + Nq + N_i ^ N. Each value of the extended index, 1 < n < (A^ + 1)^, 
thus corresponds to a unique triple. 



This includes nonphysical entries with N_i{n) < which are to be filtered out 
later. The advantage of using the extended indexing is that, as a rule, many- 
body operators written in the extended index are simple banded matrices. For 



{n - 1) mod (A^ + 1), 
n- 1- Ni{n) 



No{n) 
N-i(n) 



N + 1 
N -Ni{n)-No{n). 



(A.l) 
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example, if A; = 1, 0, 



{n' \ak\n) = 5n',n-/^k\P^k{p), = (iV + 1) 



i-fe 



(A.2) 



That is, in the extended index, oi and Oq are sparse matrices with a single 
nonzero side diagonal. Generating such matrices numerically is straightfor- 
ward; a_i which is not found in such simple way can then be obtained by a 
suitable permutation of rows and columns of oi (say). In the matrices gener- 
ated using (A.2) one should retain only columns for which A^_i(n) > and 
rows for which < N_i{n') < - 1 and Ni{n'),NQ{n') < N - I. (Recall 
that cife acts from the subspace of A^ atoms into the subspace of A^ — 1 atoms.) 
This is equivalent to building a "physical" linear index, changing A^i first, and 
A^o — second, while watching for the condition Ni + Nq < N; however the prop- 
erty of having a single nonzero side diagonal does not hold for the physical 
indexing. Generalisation to higher-order spin values is straightforward. 

There is a further subtlety in how products of field operators are implemented 
numerically. An unwanted by-product of introducing field operators is a for- 
mal non-conservation of the particle number. We wish to restrict ourselves to 
a subspace of given A^ and introduce annihilation operators as acting from 
this subspace into the subspace of A^ — 1 atoms. Then, only bilinear normally 
ordered operator products, aldi, are immediately defined for numerical pur- 
poses. Viewed literally, quartic normally ordered operators products [cf. the 
first line of (9)] imply knowledge of the field operators acting from the sub- 
space of A^ — 1 atoms into the subspace of A^ — 2 atoms and hence cannot be 
implemented under this restriction. The solution to this problem is obvious: 
use commutation relations. So, with being the unit matrix in the subspace 
of N atoms. 



etc. [cf. also the second hne of (9)]. Equations (A. 3) hold in the complete Fock 
space, but only their right-hand sides make numerical sense under the above 
restriction. 
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